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(1.1)  A  single  broadcast  channel  is  shared  among  an  infinite  number  of 
independent  sources  which  emit  packets.  Channel  time  is  discrete  and 
it  takes  one  interval  of  time  to  transmit  one  packet.  Sources  are 
synchronized  so  that  transmissions  are  initiated  at  the  beginning  of 
each  time  slot  t  =  0,1,2,...  .  We  denote  by  X  the  number  of  new  pac¬ 
kets  generated  in  a  time  slot  and  assume  that  X  has  a  Poisson  distri¬ 
bution  with  parameter  X  so  E(X)  =  X  and  variance  V(X)  =  X.  X  is 
also  independent  of  the  past  history  of  the  channel. 


(1.2)  Each  transmission  is  broadcast  to  every  user  -  including  the  emitter 
When  several  stations  transmit  simultaneously,  packets  will  collide  and 
none  of  them  is  received  correctly.  When  n  packets  collide  this  is 
called  "a  collision  of  multiplicity  n".  Each  user  involved  in  the 
collision  uses  the  same  algorithm  for  retransmitting  its  blocked  packet 
and  they  must  resolve  the  collision  without  the  benefit  of  any  other 
source  of  information  on  other  users  except  for  information  concerning  the 
channel  itself, which  provides  ternary  feedback:  idle,  successful  trans¬ 
mission,  collision. 
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The  first  algorithm  proposed  for  resolving  the  collision  is  the 
ALOHA  algorithm  which  is  now  known  to  be  unstable  and  we  refer  the  reader 
to  the  author’s  paper  {Ro-To  1983]  for  a  proof  using  elementary  martin¬ 


gale  theory.  The  first  class  of  stable  algorithms  were  developed  in¬ 
dependently  by  [Ca  1979],  [T-M  1978]  with  additional  interesting  varia¬ 
tions  due  to  Massey  [MA  1981],  and  Tsybakov-Vvedenskaya  [T-V  1980].  Of 
primary  importance  in  each  of  these  papers  is  the  determination  of 

1  =  maximum  allowable  arrival  rate  of  new  packets.  Tn  particular 

max 

for  \  <  X  the  system  is  stable  -  equivalently,  the  Markov  chain  which 

describes  the  system  is  positively  recurrent.  As  noted  by  several 

authors  including  [Ma  1981],  [FFH  1983]  the  only  quantity  one  needs  in 

order  to  evaluate  X  is  the  expected  length  of  the  collision  resolu- 
max 

tion  interval  (CRI)  determined  by  the  CRA.  More  precisely  let  H(n) 
denote  the  random  length  of  time  required  by  the  CRA  to  resolve  a  col¬ 
lision  of  multiplicity  n.  The  CTM  algorithm,  see  [Ma  1981],  leads 
to  the  recurrence  relation 


£(n)  »  1  +  2(S  )  +  H(n  -  S  ),  n  >  2 
(1.3)  n  n 

HO)  =  HU) 

where  is  a  binomially  distributed  random  variable  with  p(sn  =  j)  = 

Setting  L(n)  =  E(H(n))  and  using  the  fact  that 


E(H(n) 


n  I  _-n 


T.  E(H(n)|Sn  =  i)  [ij 2 


=  £  (1  +  L(i)  +  L(n  -  «) 

i=*0 


(")2 


1  +  2  £  L(i) 


f?l  2-n  . 


Consequently  L(n)  satisfies  the  recurrence  relation 


(1.4) 


L(n)  =1+2  Z  L(i)  "  2~n,  n  >  2 
i=0  '  > 

L(0)  =  L(l)  =  1  . 


Remark:  This  is  equation  (3.12)  on  p.  85  of  (Ma  1981]. 

Similarly  the  CTM  algorithm  with  continuous  entry  satisfies  a  re¬ 
currence  relation  of  the  form: 

£(n)  =  1  +  £(S  +  X)  +  £(n  -  S  +  Y) 

n  n 

£(0)  =  £(1)  =  1 


(1.5) 


where  S  ,  X,  Y  are  mutually  independent; 

P(X  =  j)  =  P(Y  =  j)  -  exp(-X)Xj/j! ,  j  =  0,1,... 


and  is  binomially  distributed  as  before  with  parameters  n,  p  =  1/2. 

Proceeding  in  exactly  the  same  way  as  we  did  in  deriving  (1.4)  we  see 
that  L(n)  =  E(£(n))  satisfies  the  system  of  equations  (with  infinitely 
many  unknowns ! ) 


where 

(1.7)  p  -  P(S  +  X  =  j)  =  Z  W2‘nexp(-X)Xj"1/(j-i)«;  j  >  0. 

nJ  n  i=0  I.1! 

For  future  reference  we  note  that 
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(1.8)  p  n=  2  °e  ^  and  p,=2ne^(A+n). 

nu  ni 

It  is  worth  pointing  out  that  similar  but  more  complicated  recurrence 
relations  for  the  expected  delay  have  been  derived  by  f FFHJ  1983]  and  it 
would  therefore  seem  to  be  quite  useful  to  present  an  elementary  method 
for  studying  the  existence,  uniqueness  and  asymptotic  properties  of  solu¬ 
tions  to  equations  of  the  form  (1.6).  It  is  the  purpose  of  this  paper 
to  present  just  such  a  method,  a  percursor  of  which  may  already  be  found 
in  a  previous  paper  of  the  author  [Ro  1973].  The  method  is  elementary 
in  the  sense  that  the  intricate  and  ingenious  complex  variables  methods 
of  [FFH  1983]  are  circumvented.  It  is  to  be  observed  that  a  similar 
idea  also  appears  in  [T-V  1980]  where  they  refer  to  the  notion  of  a 
"barrier  function". 

The  paper  is  organized  as  follows:  In  part  2  we  prove  an  existence 
and  uniqueness  theorem  for  equation  1.6.  In  particular  for  X  <  ^max  = 

(-5  +  AT)/4  =  .35078  we  show  that  there  exists  a  unique  nonnegative 
solution  L(n)  to  (1.6)  satisfying  the  growth  condition 

(1.9)  L(n)  <  a*n  +  b  for  n  ^  2. 

Some  additional  results  obtained  include: 

(i)  a  lower  bound  for  L(n)  of  the  form 

(1.10)  L(n)  >  [ 2 (n  -  X)  -  1]/(1  -  2X) ,  n  >  2,  X  <  X 

—  —  max 

For  example  when  X  =  .34  we  obtain 

l.(n)  >  6.25n  -  5.25,  n  >  2. 


v  ^  v  *-  - 


(1.11) 


On  the  other  hand  when  L(n)  is  the  solution  to  (1.4)  then  it 
can  be  shown,  cf.  |Ma  1981,  p.  87]  that 

(1.12)  L(n)  <  3n  for  n  >_  2. 

This  confirms  in  a  quantitative  way  one's  intuition  that  the  length 
of  a  CRI  is  longer  on  the  average  for  the  continuous  entry  scheme  than 
it  is  for  the  CRA  with  blocked  access.  The  lower  bounds  (1.10)  and  (1.11) 
appear  to  be  new. 

(ii)  A  short  and  simple  proof  of  the  intuitively  obvious  fact  than  L(n) 
is  monotone  increasing  in  n  -  see  lemma  2.19. 

(iii)  Simplifications  and  improvements  in  the  statements  and  proofs  of 
Theorem  2  of  [T-V  1980],  In  particular  their  condition  (4.5)  is  shown 
to  be  unnecessary  for  the  validity  of  their  theorems. 

In  part  3  we  turn  our  attention  to  the  sequence  L(n)  defined  by 
the  recurrence  relation  (1.4)  already  treated  in  some  detail  in  [Ma  1981]. 
In  that  paper  Massey  conjectured  (see  inequality  3.31  of  [Ma  1981])  that 

(1.13)  L(i)  +  L(n  -  i)  <  L(n)  +1,  0  <  i  <  n 

and  used  this  inequality  to  derive  upper  and  lower  bounds  for  the  variance 
V(n).  It  seems  not  to  have  been  noticed  that  the  validity  of  (1.13)  im¬ 
plies  the  existence  of  lira  L(n)/n  which  conflicts  with  assertions  in 

n-*-® 

the  literature  that  this  limit  does  not  exist.  Thus  Massey's  derivations 
of  upper  and  lower  bounds  for  V(n)  seem,  to  this  atuhor  at  least,  to 


bo  incomplete. 


2.  Existence  and  Uniqueness  of  Solutions  to  Equation  (1.6) 


To  simplify  the  tyepsetting  we  denote  a  sequence  by  f (0) ,f (1) , . . . 
instead  of  the  more  customary  notation  f^.f^,...  and  define  the 
positive  linear  operator  acting  on  sequences 

f(n),  n  =  0,1,2  via  the  recipe 

CD 

(2.1)  Gnf  £  2  ^  f(j)Pnj  =  2E(f(Sn  +  X)). 

In  this  notation  the  infinite  system  of  equations  (1.6)  can  be  re¬ 
written  in  the  more  convenient  operator  form 

L(n)  -1+GL,  n>2 

(2.2)  n 

L(0)  -  L(l)  -  1  . 

We  can  get  rid  of  the  constant  1  appearing  on  the  right  hand  side  of 

(2.2)  via  the  change  of  variables  f(n)  =  1  +  L(n),  n  j>  0  and  it  is 
easy  to  see  that  f(n)  satisfies  the  system 

f(n)  =  G  £ 

(2.3) 

f(0)  -  f(l)  -  2  . 

We  begin  by  noting  that  the  function  h(n)  *  M(n  -  2A),  n  ■  0,1,2, 
satisfies  the  equation  (2.3)  but  with  slightly  different  initial  con¬ 
ditions:  h(0)  ■  -2AM,  h(l)  -  M(1  -  2A).  To  see  this  just  compute 

G  h  -  2E(h(S  +  X))  -  2ME(S  +  X)  -  4MA 
n  n  n 

=  2M((n/2)  +  A,  -  4MA  -  M(n  -  2A)  -  h(n). 

Note:  Throughout  this  paper  we  assume  0  <  A  <  1/2.  Notice  that  h(n) 
is  linear  in  n  and  it  is  reasonable  to  conjecture  that  f(n)  itself 
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is  nearly  linear  in  the  sense  that  f (n)  a*n  +  b  as  n  -*■  Unlike 
tlie  system  of  equations  (1.4)  for  which  uniqueness  anti  existence  is 
trivial  it  is  not  at  all  a  priori  evident  that  solutions  to  (2.3)  are 
unique.  We  begin  therefore  with  a  proof  of  uniqueness  since  this  re¬ 
quires  the  introduction  of  a  technique  that  will  serve  us  well  in  our 
existence  proof. 

Theorem:  There  is  at  most  one  nonnegative  solution  f(n)  to  (2.3) 
satisfying  the  growth  condition  |f(n)|  <_  a*n  +  b. 


Proof:  Suppose  f^(n)  and  f2(n)  are  two  solutions  to  (2.3)  then 
W(n)  =  f^n)  -  f2(n)  satisfies 


(2.4) 


w(n)  =  G  w,  n  >  2 
n  — 

w(0)  •  w(l)  ■  0  and  |w(n) |  £  a'n  +  b' 


Thus  it  suffices  to  show  that  the  only  solution  to  (2.4)  is  w(n)  =  0. 

This  follows  at  once  from  the  existence  of  a  "barrier  function" 

2 

p(n)  =  h(n)  +  n  ,  n  >  2,  p(0)  =  p(l)  “  0  with  the  following  pro¬ 

perties: 


(i)  p(n)  >0,  n  >  2 

(2.5)  (ii)  lim  |w(n)|/p(n)  -  0 

n-*  00 

(iii)  2  Ep(j)p  .  *  ITT',"*  p(n)  n  >  2  . 

j»2  n 


-2 

(i)  is  obvious  and  (ii)  follows  from  the  fact  that  lim  n  w(n)  *  0, 

n+oo 

while  1  -  2A  >  0.  Later  we  shall  have  to  restrict  A  even  further 


X 


and  assume  \  ■»  \ 


=  . 35078.  We  postpone  for  the  moment  a  proof  of 


max 

(iii)  and  proceed  to  derive  the  uniqueness  theorem  by  exploiting  the 
properties  of  p.  Lndeed  the  assumption  that  w(n)  t  0  leads  to  a 
contradiction  as  we  now  show.  Condition  (ii)  implies  that  if  w(n) 
is  not  identically  zero  then  there  exists  a  finite  integer  n'  2 

with  the  property  sup  jw(n)|/p(n)  =  | w(n' ) | /p(n ' )  0.  From  equation 

n>2 

(2.4)  we  see  that 

| w (n ’ ) |  <  2  I  l |w(j) |/p(j)]p(j)p  , 
j  =  2  n  3 

<  (iw(n') |/p(n'))Gn,p  <  |w(n')| 

where  we  have  used  (2.5  Iii)  in  the  last  step.  This  is  a  contradiction 

so  we  conclude  sup  |w(n)|/p(n)  =  0  i.e.  w(n)  =  0. 
n>2 

We  now  prove  that  p(n)  satisfies  property  (2.5  iii).  Since 

oo 

..  2 

2  1.  i  p  .  is  twice  the  second  moment  of  the  random  variable  S  +  X 

j  =  l  nj 

03  2 

whose  variance  is  equal  to  n/4  +  \  we  see  at  once  that  2  ?.  j  p  ,  = 

2  nj 
n/4  +  A  +  (n/2  +  A)  .  Consequently 

00  2  2 

(2.6)  2  Y  Jp  =  2[  (n/4  +  A)  +  (n/2  +  A)Z]  -  2pnl> 

j=2  J 

pnl  -  2"VA(X  +  n),  cf  (1.8). 

A  similar  calculation  yields 

(2.7)  2  *  h(j)p  -  h(n)  -  2p„M(-2A)  -  2P  ,M(1  -  2A)  . 

j_2  nJ  nU  ni 

After  some  routine  algebraic  manipulations  the  term  4MXp  n  -  2M(1  - 

nu 

simplifies  to  M2  n+^exp (- X) ( 2X ^  +  X(2n  +  1)  -  n).  A  necessary  condition 


2X)P) 


then  for  (2.5  iii)  to  hold  is  that  2A  +  A(2n  +1)  -  n  ^  0  for  all 

n  2.  This  requires  an  examination  of  the  largest  positive  root 

2 

A^ax^n)  °f  the  quadratic  equation  2A  +  \(2  n  +  1)  -  n  which  turns 
out  to  be 


(2.8) 


A  (n)  =  f  —  ( 2n  +  1) 
max 


+  1)  +  An2  + 


12n  +  1  ]/4  . 


2  2 

Notice  that  4n  +  12n  +  1  >  (2n  +  2)  for  n  _>  2  and  therefore 

•\  (n)  >  .25.  Actually  we  can  do  a  lot  better  since  it  is  easy  to  see 

that  the  smallest  value  of  A  (n)  occurs  at  n  =  2  with  A  (2)  = 

max  max 

.35078  =  A  .  Summing  up  then  we  have  shown  that  for  A  <  .35078 
max 


(2.9a)  2  I  h(j)p 

J-2 


nj 


h(n)  <  M2_n+1exp(-A)(2A2  +  X(2n  +  1)  -  n)  <  0. 


The  next  step  is  to  study 
we  see  that 


2  £  j2p  -  n2]  for  n  >  2.  From  (2.6) 


j-2 


“J 


(2.9b) 


2  I  j2p  -  n2  =  (n/2  +  2A  +  n2/2  +  2nA  +  2A2  -  n2 
J-2  nj 

=  n/2  +  2A  +  2nA  +  2A2  -  n2/2. 


A  simple  calculation  shows  that  the  right  hand  side  of  (2.10)  is  strictly 
negative  for  all  n  >  4  and  0  £  A  <.  .5.  And  it  is  also  clear  that 
for  any  fixed  A  <  Aq  and  n  -  2,3  we  can  choose  M  (depending  on  A) 

CO 

sufficiently  large  and  positive  so  that  2  L  (h(j)  +  j2)p  <  h(n)  +  n2,n  >  2 


j-2 


nj 


and  this  completes  the  proof  of  uniqueness. 


Remark:  The  function  p(n)  -  h(n)  +  n  for  n  s  2  with  p(0)  -  p(l)  -  0 


satisfies  the  inequality  G  p  <  p(n)  and  is  called  a  "barrier  function" 


by  Tsybakov-Vvedenskaya .  By  constructing  barrier  function  P(n)  ten- 

2 

ding  to  infinity  at  the  rate  n  the  uniqueness  proof  of  Tsybakov- 
Vvedenskaya  is  considerably  simplified  as  the  reader  may  easily  check 
for  himself.  Similar  ideas  were  used  (but  with  respect  to  a  different 
operator  G!)  in  our  paper  [Ro  1973]. 

We  now  prove  the  existence  of  a  nonnegative  solution  to  the  system 
of  equations 


(2.10) 


Gng  =  g(n) 

g(0)  = aQ,  g(l)  =  a a.  0,  i  *  0,1. 


Remark:  The  system  (2.3)  corresponds  to  the  special  case  a^  =  a^  =  2. 
Our  proof  is  similar  to  that  in  [T-V  1980]  except  that  a  judicious  use 
of  the  Lebesgue  dominated  convergence  theorem  yields  a  more  streamlined 
proof.  We  need  the  following 


(2.11)  Lemma:  Let  X°(n)  =  h(n)  ,  n  2,  and  X°(0)  =  a^,  X°(l)  =  a 
Then  for  A  <  Aq  X°(n)  is  a  barrier  function  i.e.  Gnx°  £  X°(n) , 
n  _>  2,  provided  M  =  M(A,aQ,a^)  defined  at  (2.15)  below. 


Proof:  The  necessary  calculations  are  similar  to  those  already  carried 
out  at  (2.7).  Since  X°(j)  -  h(j)  ^  0  only  for  j  =  0,1  we  see  that 

Gn(X°  -  h)  =  (X°(0)  -  h(0))pnQ  +  (X°(l)  -  h(l))pnl  =  r(n.V). 


Thus  for  n  >  2  we  have 


(2.12)  G  X°  =  G_h  +  G_(X°  -  h)  «  X°(n)  +  r(n,A),  n 


Consequently  r(n,A)  ^  0,  n  £  2  is  a  necessary  and  sufficient  condition 
for  X°  to  be  a  barrier  function.  In  particular 

(2.13)  r (n, A)  =  M2_ne“X(2A2  +  A(2n  +  1)  -  n)  +  aQ  e~A2~n 

+  a^e  A\2  n(A  +  n) 

=  2_ne‘A{M(2A2  +  \(2n  +  1)  -  n)  +  aQ\  +  a^A(A  +  n]  . 

Since  a^  and  a^  are  nonnegative  it  is  clear  that  a  necessary  con 

2 

dition  for  r(n,A)  <_  0  all  n  2  is  that  2A  +  A(2n  +  1)  -  n  <  0  - 

equivalently  for  A  <  A~  =  A  (2)  =  .35078  (cf  (2.8)).  Keeping  in 

u  max 

2 

mind  that  n  -  (2A  +  A(2n  +1))  >0  for  A  <  A^  we  choose 

(2.14)  M  >  (a  An  +  aQ  +  axA2)/(n  -  (2A2  +  A(2n  +  1))),  n  _>  2, 

and  it  is  easy  to  see  that  the  right  hand  side  of  (2.14)  remains  bounded 
as  n  -*■  °°.  So  it  suffices  to  choose 

(2.15)  M(Apa_,a  )  =  sup  [(a.An  +  a_A  +  a1A2)/(n  -  ( 2 A 2  +  A(2n  +  1)))J. 

1  n>2  1 

This  completes  the  proof  of  the  existence  of  a  barrier  function. 

Remark:  Notice  that  Gnx°  i  X°(n)  implies  that  X°  is  a  nonnegative 

function  integrable  in  the  sense  of  Lebesgue  with  respect  to  the 

GO 

"weights"  bn(j)  =  2p  .  since  G^X0  *  A  X°(j)p  (j)  £  X°(n).  We  shall 
o  nj  n  j=0 

make  use  of  this  remark  shortly. 

Equipped  with  the  barrier  function  X°(n)  of  lemma  (2.11)  it  is 
fairly  straightforward  to  construct  a  solution  to  equation  (2.10).  De- 
fine  the  sequence  of  functions  X  (n)  via  the  recipe 


6 


Xk(0)  =  a 

(2.17)  Assertion:  0  £  X 

lim  X  (n)  =  g(n)  exists, 
k  <■ ->■ 


o-  xk<>>  ■ 


k+1 


(n) 


Xk(n) 


Moreover  G 

n 


,  k  =  0,1,2,. 

g=  g(n)  for 


n 


and  therefore 
>  2  and 


g(n)  <.  X°(n)  =  M(X,aQ,a1)(n  -  2X)  ,  so  jg(n>!  <_  a*n  +  b. 

['roof:  X*(n)  =  GnX0  <  X°(n)  implies  X2(n)  =  G  X1  <  G^0  =  X1  (n)  <  X°(i 

k 

and  so  on.  Thus  X  (n)  is  monotone  decreasing  in  k  as  claimed  and 

Ic 

since  the  sequence  is  bounded  below  it  follows  that  ft(n)  *  lim  X  (n) 

k+oo 


exists  and  satisfies  the  inequality  |g(n)|  <_  a*n  +  b.  The  only  thing 
we  have  to  check  is  that  G^g  =  g(n)  for  n  _>  2.  Observe  that 

,  OO 

0  -  X  (j)  X°(j)  and  I  X°(j)2p  £  X°(n)  permit  us  to  use  the 

j=0  nJ 

Lebesgue  dominated  convergence  theorem  to  conclude 


(2.18)  lim  X  Xk(j)2p  -  l  (lim  Xk(j))2p 
k  +  lJO  j=0  J  j=0  k-*- °°  J 

i.e.  lim  GnXk  =  G  (limXk)  =  G  g.  But  |g(n)  -  G  g |  <_  |g(n)  -  Xk+1(n)| 

|xk+1(n)  -  Gng|. 

Clearly  lim  |g(n)  -  Xk+^(n)|  *  0  and  lim  |Xk+*(n)  -  G  g|  = 

co  k-*-  00 

lim  | G  Xk  —  G  g |  =0  using  (2.18)  for  the  last  step.  This  completes 

k  >•« 

the  proof  of  the  existence  of  a  solution  to  (2.10)  satisfying  the  growth 
condition  |g(n)|  a*n  +  b,  consequently  the  solution  is  unique. 

We  now  proceed  to  use  the  existence  and  uniqueness  theorems  to  de¬ 


rive  the  lower  bound  (1.10). 


Set  w(n)  =  f(n)  -  h(n)  where  f  Is  the  unique  solution  to  (2.3) 
satisfying  the  growth  condition  ] f (n) |  <  an  +  b  and  h(n)  =  M(n  -  2X). 

It  is  easy  to  check  that  for  0  M  <■  2(1  -  2X)  *  we  have  w(0)  >  0, 
w(l)  >0  with  G^w  =  w(n),  n  2  and  |w(n) t  £  an  +  b.  Consequently 

by  the  uniqueness  theorem  w(n)  >  0  all  n  >  2  i.e. 

f (n)  >  h(n)  =  [2(1  -  2X)-1  -  e](n  -  2X)  for  n  >  2. 

Since  this  is  true  for  every  r  >  0  we  get  f(n)  >  (2(1  -  2X)  *)(n  -  2X). 

The  lower  bound  for  L(n)  is  now  derived  by  noting  that  L(n)  «  f(n)  -  1. 

Finally  we  show  how  the  probabilistic  representation  GRg  =  2E(g(Sn  +  X)) 
may  be  used  to  show  that  the  solutions  to  (2.10)  are  monotone  increasing 
in  n. 

(2.19)  Lemma:  If  g(n)  is  monotone  increasing  and  nonnegattve  then  so 
is  Gng. 

Proof :  Cn+1g  -  Gng  =  2E(g(Sn+1  +  X)  -  g(Sp  +  X)},  so  it  suffices  to 
show  that  E(g(Sn+1  +  X)  -  g(Sn  +  X))  >  0.  But 

E(g(Sn+1  +  X)  -  g(Sn  +  X)  i  S^  =  i  ,X  =  .1)  = 

(1/2  ) (g(i  +  j  +  1)  -  g(i  +  j))  >_  0  for  every  i,j. 

Equivalently 

(2.20)  E(8(Sn+l  +  X>  '  8(Sn  +  X> | Fn)  >  0 

where  F  =  smallest  sigma-field  with  respect  to  which  both  S  ,X  are  measurable 
n  n 

Taking  expectations  of  both  sides  of  (2.20)  yields  the  desired  inequality 

k 

that  ^’n+^R  "  1.  0.  In  particular  then  each  function  X  (n)  is  mono¬ 

tone  increasing  in  n  since  X°(n)  is  obviously  monotone.  Consequently 
so  is  the  limit  g(n)  =  lim  X^(n) . 


J.  Remarks  on  the  asymptotic  behaviour  of  the  solution  to  (1.4) 

In  the  postscript  to  iMa  1981]  reports  that  Vvedenskaya  has  shown 

(so  far  unpublished)  that  lim  l.(n)/n  does  not  exist  although  it  can  be 

n->-  a> 

shown  (and  we  shall  sketch  such  a  proof  below  based  on  the  methods  in 
the  preceding  section  of  this  paper)  that 

(3.1)  2.881  <  lim  inf  L(n)/n  <  lim  sup  L(n)/n  <  2.8966. 

n-*  00  n-* on 

In  the  same  paper  Massey  asserted,  without  proof,  that 

(3.2)  L(i)  +  L(n  -  i)  <  L(n)  +1,  0  <  i  <  n. 

In  September  of  1983  I  pointed  out  to  my  colleagues  Don  Towsley  and 
.Jack  Wolf  of  the  Department  of  Electrical  and  Computer  Engineering  here 
at  the  University  of  Massachusetts  that  any  inequality  of  the  form: 

(3.3)  L(i)  +  L(n  -  i)  <  L(n)  +  c,  0  <  i  <  n  where  c 

is  an  arbitrary  constant,  necessarily  implies  the  existence  of  lim  L(n)/ 

n-*.  ™ 

In  other  words  the  nonexistence  of  the  lim  L(n)/n  is  inconsistent 

n-*'  '*• 

with  the  validity  of  (3.3).  To  see  this  we  need  to  introduce  the  notion 
of  a  subadditive  sequence. 

(3.4)  Definition:  We  say  that  a  sequence  a^  is  subadditive  if 

a  +  a  >  a  .  holds, 
n  m  —  n+m 

(3.5)  Theorem:  If  a  is  a  subadditive  sequence  then  lim  a  /n  exists 

n  n 


Proof:  See  [P-SZ  (1972)  p.23,  problem  98]. 
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Application:  Suppose  (3.1)  holds  then  the  sequence  a  *  c  -  L(n)  satis¬ 

fies  the  condition  a  +  a  >  a  which  is  equivalent  to  the  condition 

m  n-m  --  n 

a  +  a  >  a  .  i.e.  a  is  subadditive  and  therefore  lira  a  /n  =  lira  L(n)/n 
m  n  —  m+n  n  n 

n  +  ■»  n  ►  ■ 

exists. 

Now  inequality  (3.2)  was  used  by  Massey  to  derive  upper  and  lower 
bounds  on  the  asymptotic  behaviour  of  the  variance  of  the  CRI  length 
V(n)  defined  by  (1.3).  These  proofs  must  now  be  regarded  as  incomplete. 

Turning  now  to  the  proof  of  (3.1)  we  shall  content  ourselves  with  a 
sketch  of  the  derivation  of  the  upper  bound 


(3.6)  lim  sup  l.(n)/n  <  2.8966. 

n->  w 


As  we  did  in  part  2  of  this  paper  we  set  f(n)  =  L(n)  +  1 
the  asymptotic  behaviour  of  the  recurrence  relation 

f 


(3.7) 


f(n)  =  2E(f(Sn)),  n>2 

f(0)  =  f (1)  =  2  . 


and  study 


Wo  shall  make  extensive  use  of  the  following  elementary  lemma  the  proof 
of  which  is  left  to  the  reader. 


(3.8)  Lemma:  Suppose  W(i)  >0,  0  <  i  <  m  and  W(n)  >  2E(W(Sn)),  n  ^  m. 

Then  W(n)  >  0  all  n. 

Bring  in  the  function  H(n)  =  a.n  which  satisfies  the  recurrence 

relation  (3.7)  for  all  n.  Set  g(n)  =  H(n),  n  >  4  and  g(n)  =  f(n), 

n  <  3  and  compute  E(g(S  )  -  H(S  ))  =  r(n).  Using  the  fact  that  f ( 0)  = 

n  n 

f ( 1 )  =  2,  f ( 2 )  3  6  we  see  that 


=  2  "n(n  -  l)  2/ (n(n  -  1))  +  (2  -  a)/(n  -  1)  +  (3  -  0 
+  (n  -  2) (13/9  -  a/2)  . 


The 

problem 

then  is  to  choose  a 

in  such  a  way  that  r(n)  <  0 

for  n 

and 

W(n)  = 

g(n)  -  f(n)  >  0  for 

0  <  n  <  4.  From  (3.9)  we  see 

that  a 

necessary  condition  is  for  13/9  -  a/2  <0  or  a  >  2. 8889.  It  turns  out 
■i  =  2.8889  is  too  small  and  that  we  must  choose  a  =  23/9  +  .00778  =  2.8966 
With  this  choice  of  a  one  can  show  by  means  of  a  straightforward  but 
tedious  calculation  that  W(n)  =  g(n)  -  f(n)  satisfies  the  conditions 
of  lemma  3.8  and  hence  f (n)  <_  g(n)  =  2.89666n  for  n  >  4.  In  a  similar 
calculation  we  can  derive  Massey's  lower  bounds  (Table  3.2  on  P.  87  of 
I Ma  1981]). 


4.  Concluding  remarks. 


In  this  paper  we  have  presented  an  operator  method  for  obtaining 
upper  and  lower  bounds  for  the  expected  length  L(n)  of  a  CR1  for  various 
protocols.  Although  we  shall  not  carry  out  such  an  analysis  here  it  is 
noteworthy  that  the  same  methods  can  be  used  to  obtain  upper  and  lower 
bounds  for  the  expected  delay  W(n)  =  expected  delay  experienced  by  a 
packet  that  had  its  first  transmission  during  a  collision  of  multiplicity 
n.  As  pointed  out  by  [FFHJ  (1983)]  in  the  case  of  CTMCRA  with  blocked 
access  w(n)  (with  W(n)  =  E(w(n)))  satisfies  the  recurrence  relation 


(4.1) 


w(n)  =  w(S  .  +  1)  with  probabilitv  1/2 
n- 1 

w(n)  *  w(n  -  Sn_j)  +  ?(Sn_1)  with  probability  1/2. 


Taking  expectations  as  before  we  are  led  to  a  more  complicated  svsLem  of 
recurrence  equations  whose  asymptotics  can  be  studied  by  means  of  similar 
methods.  We  intend  to  pursue  these  questions  in  another  paper. 

We  must  however  point  out  that  our  methods  have  so  far  been  unable  to 
give  a  proof  of  Massey's  upper  and  lower  bounds  for  V(n)  =  Variance  (£(n)). 
Obtaining  upper  and  lower  bounds  for  the  second  and  higher  order  moments 
of  £(n)  is  one  of  the  more  interesting  and  still  open  problems  in  this 
area. 
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